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Abstract-  A  compilation  of  54  sets  of  geoacoustic  measurements  of  marine  sediments  was  reviewed  in  order  to  bound  estimates  of 
parameters  for  use  in  Buckingham's  grain  shearing  (GS)  theory  of  acoustic  propagation.  These  data,  for  unconsolidated  sands 
(siliciclastic  and  carbonate)  with  grain  sizes  between  0.0156  mm  and  0.57  mm,  were  all  made  in  shallow  water  (maximum  depth  of  60  m) 
sites  at  diverse  locales.  In  each  data  set,  measurements  of  the  speed  and  attenuation  of  the  compressional  wave  at  high  frequency,  the 
shear  wave  speed  at  1  kHz,  porosity,  and  bulk  density  are  sufficient  to  calculate  the  three  free  parameters  in  GS  theory  (a  material 
exponent,  a  compressional  coefficient  and  a  shear  coefficient).  The  spread  of  the  values  calculated  for  GS  parameters,  combined  with  a 
sensitivity  analysis,  do  not  support  the  use  of  a  single  material  exponent  value  for  all  sediments.  Nor  do  they  support  the  notion  that  the 
material  exponent  be  1,  which  would  be  the  case  if  elastic  and  viscous  forces  at  grain  contacts  are  equal.  Finally,  it  is  suggested  that, 
given  the  sensitivity  to  spot  measurements,  the  GS  parameters  be  estimated  simultaneously  from  a  full  range  of  compressional  and 
shear  wave  measurements. 


I.  Introduction 

Buckingham’s  Grain  Shearing  (GS)  theory  was  introduced[l]  as  an  alternative  to  poro-elastic  theories,  such  as  Biot’s[2],  in  an 
attempt  to  better  match  observed  dispersion  curves.  The  foundation  of  this  alternative  theory,  that  rigidity  in  the  sediment  is 
provided  by  grain- to-grain  shearing,  removes  the  necessity  of  describing  an  elastic  frame,  an  artificial  construct  which  must  be 
described  by  unobservable  elastic  parameters.  However,  GS  theory  introduces  three  parameters,  which  are  not  in  themselves 
observable,  but  are  derived  from  other  physical  observations.  Buckingham  presented  an  expanded  theory,  designated  the  viscous 
Grain  Shearing  theory  (VGS)  [3],  to  account  for  the  effects  of  viscosity  in  the  pore  fluid  that  affects  dispersion  curves  at  low 
frequencies. 

GS  theory  describes  stick-slip  events  at  grain  contacts,  which,  rather  than  being  truly  impulsive,  have,  on  average,  some 
exponential  decay,  parameterized  by  a  material  exponent  n,  Compressional  and  shear  waves  propagate  via  these  grain  contacts 
with  relative  intensity  parameterized  by  two  elastic  moduli,  Yp  and  Ys  respectively.  These  are  not  elastic  moduli  relating  stress  and 
strain  rates,  but  are  related  to  the  rate  and  intensity  of  radial  and  translational  shear  events  at  grain  contacts.  Not  directly 
measureable,  they  are  derived  from  observed  quantities. 

Given  a  known  porosity  A,  GS  theory  requires  estimation  of  three  unknown  parameters  {n  ,  Yp  Vs)-  Buckingham  has 
described  a  method  of  calculating  these  three  values  from  measurements  of  compressional  and  shear  wave  speeds  and 
attenuations[3].  First  compressional  sound  speed  and  attenuation  at  a  high  frequency  is  used  to  calculate  the  material  exponenent 
n.  Then  the  shear  speed  at  a  low  frequency  is  used  to  calculate  both  moduli,  Yp  Ys-  This  method  avoids  use  of  the  shear  wave 
attenuation,  a  usually  imprecise  measurement  seldom  available.  VGS  theory  adds  a  fourth  parameter,  a  viscoelastic  time  constant, 
which  is  constant  for  all  sediments.  However,  inclusion  of  this  parameter  slightly  changes  estimates  of  the  three  GS  parameters. 

II.  Geoacoustic  data  set 

The  geoacoustic  measurements  were  all  made  in  situ  using  various  versions  of  the  In  Situ  Sediment  geoAcoustic 
Measurement  System  (ISSAMS).  The  measurement  system  is  described  in  Chapter  5  of  Jackson  and  Richardson  [4]  and  most  of 
the  geoacoustic  and  physical  property  data  used  in  these  analyses  can  also  be  found  in  that  reference.  The  data  base  can  be  traced 
back  to  the  original  sources  using  references  in  [4].  All  in  situ  geoacoustic  measurements  were  made  at  20-30  cm  below  the 
sediment-water  interface  over  path  lengths  ranging  from  30-100  cm.  Compressional  wave  speed  and  attenuation  were  measure  at 
either  38  or  58  kHz  using  time-of  flight  and  amplitude  of  a  5  to  10-cycle  pulsed  sine  waves  propagating  between  identical 
radially-poled  ceramic  cylinders  through  sediment  and  a  reference  of  seawater  just  about  the  sediment  water  interface.  Shear 
wave  speed  was  measured  at  1  kHz  using  time-of- flight  between  bimorph  ceramic  benders  mounted  on  the  same  diver-deployed 
or  remotely-operated  hydraulic  systems.  Shear  wave  attenuation  was  measured  at  selected  sites  using  a  4-transducer  transposition 
technique  which  calculates  shear  wave  attenuation  from  waveform  amplitudes  measured  using  two  transmit  and  two  receiver 


0-933957-38-1  ©2009  MTS 


transducers.  This  technique  eliminates  the  need  to  measure  transducer  sensitivity  or  measure  variable  insertion  losses.  Multiple 
sediment  cores  were  collected  from  each  site  to  measure  high  frequency  (400  kHz)  compressional  wave  speed  and  attenuation 
and  to  provide  data  of  sediment  physical  properties  such  as  grain  size  distribution,  sediment  bulk  density  and  porosity.  All  values 
of  wave  speed  and  attenuation  as  well  as  bulk  sediment  physical  properties  used  in  these  analyses  are  averages  from  multiple 
deployments  and  multiple  sediment  cores  collected  at  the  same  location  (a  roughly  25-m^  area). 

The  variability  among  measured  wave  speeds  and  attenuation  at  a  single  location  is  generally  thought  to  be  equal  or  greater 
than  the  actual  measurement  error  which  is  less  than  1%  for  wave  speeds  and  less  than  10%  for  attenuation. 

For  this  study  we  restrict  ourselves  to  data  from  the  54  sandy  sites  where  the  microscopic  stress  relaxation  mechanisms  that  are 
part  of  the  Buckingham  theory  are  most  likely  to  be  applicable.  It  seems  imlikely  that  this  particle-to-particle  stress  relaxation 
mechanism  can  be  applied  to  muddy  sediment  where  electrostatic  repulsive  and  attractive  forces  and  the  adhesion  of  organic 
matter  control  particle-to-particle  interactions  and  the  flexure  of  clay  particles  may  provide  a  dissipation  mechanism. 

III.  SENSITIVITY  ANALYSIS 

Errors  in  estimates  of  Buckingham’s  parameters  are  driven  by  uncertainty  in  the  measurements  of  the  physical  parameters  from 
which  they  are  calculated.  Given  the  high  correlation  between  porosity  and  grain  size,  only  one  need  be  included,  and  here  it  is 
porosity.  Buckingham  describes  a  process  of  calculating  all  parameters  by  first  determining  co  as  a  function  of  the  measured 
porosity  from  Wood’s  equation[5],  calculating  the  spectral  exponent  n  from  and  Op^  (measured  at  a  high  frequency,  either 
38  or  58  kHz)  and  finally  evaluating  the  GS  moduli  from  the  shear  wave  speed  c,'”  (measured  at  a  low  frequency,  here,  1  kHz). 
Hence,  the  error  in  estimating  n  involves  errors  in  Cp"*,  and  Estimates  of  the  GS  moduli  additionally  involve  errors  in 
measuring  Assuming  the  errors  in  measurements  are  Gaussian  and  uncorrelated,  the  overall  error  in  each  estimate  is  the  sum 
of  the  contributions  from  individual  sources. 
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Here  Oq  is  the  variance  of  the  estimate  of  the  parameter  0.  The  partial  derivatives  given  in  equation  (1)  are  given  in  the  annex. 
Although  the  variance  of  the  errors  in  the  measurement  of  physical  and  geoacoustic  parameters  vary  for  each  data  set,  for  this 
analysis,  it  was  assumed  that  the  coefficient  of  variation  in  each  measurement  was  the  same  for  all  data  sets.  Hence  the  following 
variances  were  assumed. 
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Errors  in  estimates  of  the  GS  model  in  turn  result  in  errors  in  the  dispersion  curves  upon  which  they  are  based.  Only  the  errors 
for  the  compressional  wave  parameters  are  given  below.  Similar  expressions  for  shear  wave  parameters  are  trivially  similar. 
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IV.  RESULTS 

The  values  calculated  for  the  material  exponent  a7,  compressional  modulus  and  shear  modulus  ^  are  shown  in  Figs.  1  to  3, 
respectively.  Fifteen  data  sets  yielded  negative  values  for  an  unrealistic,  yet  mathematically  possible  result,  and  are  hereafter 
omitted  from  the  analysis.  For  each  value  calculated  from  the  remaining  39  data  sets,  error  bars  span  plus  or  minus  one  standard 
deviation  about  the  estimate.  Each  parameter  is  plotted  as  a  function  of  porosity.  Carbonate  sands  are  plotted  in  cyan, 
siliciclastic  in  blue,  and  individual  SAX99  data  in  red.  Values  derived  from  SAX99  used  by  Buckingham  are  green[6]. 


Porosity 


Figure  1  -  Material  exponent  versus  porosity. 


Although  there  is  a  wide  spread  in  uncertainty  in  the  material  exponent  shown  in  Fig.  1,  there  appears  to  be  a  dependence  on 
porosity.  The  bulk  of  the  values  are  of  magnitude  much  less  than  one.  This  contradicts  the  notion  that  elastic  and  viscous  forces 
at  the  grain  contacts  should  be  considered  of  equal  importance,  at  least  within  context  of  the  GS  theory.  The  apparent  increase  in 
value  for  n  with  increasing  porosity  contradicts  Buckingham's  assertion  that  one  value  for  w  can  be  used  for  all  sediments[7]. 


Porosity  is  the  determining  factor  in  evaluating  the  low  frequency  limit  of  compressional  sound  speed  c©,  given  by  the  Wood 
equation.  Co,  in  turn,  is  used  in  calculating  n.  But  there  is  obviously  further  dependence  of  n  on  porosity. 

Assuming  that  estimates  of  porosity  and  compressional  sound  speed  are  accurate,  and  that  inaccuracies  in  measuring 
compressional  attenuation  prevent  accurate  determination  of  «,  an  attempt  was  made  to  find  a  single  value  of  n  that  resulted  in 
values  of  Op  that  is  consistent  with  the  entire  data  set.  Specifically,  if  is  the  attenuation  required  to  give  a  hypothetical 
material  exponent  and  a^!  is  the  measured  attenuation  from  data  set  /,  then  a  residual  error  can  be  computed. 

^=a‘p-(/p 

(4) 

The  value  of  r/'  that  minimized  the  mean  square  error  for  all  data  sets  is  n'  =.14.  Adopting  a  lower  value  of  n'  for  all  sediments 
{n^  =.0851),  as  Buckingham  has  suggested[3],  implies  that  compressional  attenuation  used  in  the  GS  model  is  significantly  lower 
than  that  observed  in  these  data  sets.  However,  this  is  consistent  with  his  assertion  that  GS  theory  tends  to  giver  lower  bounds  on 
attenuation  values,  in  that  it  accounts  only  for  intrinsic  attenuation[3].  Other  sources,  such  as  scattering  of  sound  by  large-scale 
inhomogeneities  are  not  accounted  for  in  this  theory. 

Both  grain-shear  moduli,  plotted  in  Figs.  2  and  3,  also  show  strong  dependence  on  porosity.  Note  the  errors  are  relatively 
greater  for  carbonate  sands  than  si lici clastic.  Also  the  error  in  the  estimate  of  is  higher  than  that  for  yp.  Generally,  the  largest 
source  of  error  in  the  estimate  of  Yp  is  due  to  errors  in  measurement  of  while  that  for  the  estimate  of  ^  is  due  to  errors  in 
measurement  of 

The  uncertainty  in  estimates  of  the  GS  parameters  carries  through  to  estimates  of  sound  speed  and  attenuation  (both 
compressional  and  shear).  For  instance,  the  GS  predictions  of  velocity  ratio  (compressional  speed  to  speed  of  sound  in  water)  and 
attenuation  for  the  SAX99  data  set  used  in  [6]  are  shown  in  Figs.  4  and  5,  respectively.  The  GS  estimates  are  shown  in  blue  while 
the  upper  and  lower  bounds  (+/-  one  standard  deviation),  as  given  by  equation  (3),  are  shown  in  green  (the  lower  bound  of  the 
attenuation  is  zero  for  all  frequencies).  Also  shown  for  reference  are  the  VGS  estimates  in  red.  The  error  bounds  are  based  on 
partial  derivatives  of  the  GS  theory.  It  is  expected  that  similar  partial  derivatives  for  the  VGS  theory,  although  much  more 
complicated ,  will  differ  little  from  those  for  the  GS  theory. 


Figure  4  -  Compressional  wave  speed  predicted  by  GS  (blue)  and  VGS  (red)  theories  using  SAX99  measurements.  Error 
in  GS  theory  bounded  by  green  lines. 
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Figure  5  -  Compressional  wave  attenuation  predicted  by  GS  (blue)  and  VGS  (red)  theories  using  SAX99  measurements. 
Error  in  GS  theory  bounded  above  by  green  line,  below  by  0. 


The  uncertainties  given  in  the  above  figures  are  calculated  for  the  method  of  evaluating  the  three  GS  parameters  specified  by 
Buckingham.  However,  it  is  postulated  that  alternative  methods  may  yield  more  accurate  estimates.  For  instance,  rather  than 
sequentially  determining  n,  then  the  GS  moduli,  they  may  all  be  fit  to  observations  simultaneously.  The  possible  reduction  in 
uncertainties  in  these  estimates  is  discussed  next. 


V.  BOUNDS  ON  ALTERNATIVE  ESTIMATES  OF  BUCKINGHAM’S  PARAMETERS 


In  a  mapping  G  of  a  set  of  N„  model  parameters  to  a  set  of  observations  d, 
d  =  G(m) 

(5) 

the  covariance  of  the  A  Posteriori  errors,  C’m,  in  the  estimates  of  m  can  be  calculated  if  all  errors  are  assumed  to  be  Gaussian 
and  the  mapping  can  be  linearized  about  some  estimate  mo.  In  this  case [8], 


C'„=(c„+G[(CT  +  CDr'G,)'' 


(6) 


Here  C’^  is  the  covariance  of  the  A  Priori  errors,  Cp  is  the  covariance  of  the  data  (measurement  errors),  Ct  is  the  covariance  of 
theoretical  errors  (due  to  mismatch  between  the  forward  model  and  reality),  and  Gi  is  given  by  the  linearization 


In  this  case,  m  is  the  set  of  GS  parameters  (w,  Yp,  d  is  the  set  of  observations  (A^,  Cp,  Cs),  the  A  Priori  errors  are  given  by 
the  previous  sensitivity  analysis  and  the  linearization  provided  by  the  partial  derivatives  previously  derived.  Measurement  errors 
are  assumed  to  be  dominant,  and  theoretical  errors  are  neglected. 

In  the  case  in  which  the  three  GS  parameters  are  estimated  simultaneously,  the  standard  deviation  of  the  error  in  estimating  n 
(for  the  composite  SAX99  data  set)  is  only  slightly  reduced  (by  a  factor  of  .95),  while  those  for  the  two  GS  moduli  are  roughly 
halved.  This  is  understandable,  as  Buckingham's  sequential  method  uses  the  most  reliable  information  to  first  calculate  n,  and  the 
GS  moduli  are  based  on  a  more  error  prone  measurement  of  the  shear  wave  speed. 

Another  alternative  was  investigated,  to  see  if  uncertainties  can  be  reduced  more  significantly.  A  second  set  of  compressional 
wave  measurements  from  the  SAX99  experiment,  provided  by  the  diver  deployed  “attenuation  array"  gives  Cp  and  at  100 
kHz[9].  Assuming  similar  variances  in  these  estimates  as  with  the  58  kHz  data,  the  above  error  analysis  again  resulted  in  an 
insignificant  decrease  in  the  error  in  estimating  w,  about  a  halving  of  the  error  in  estimating  but  a  reduction  in  the  error  in  Yp  by 
a  factor  of  .07. 


VI.  CONCLUSIONS 

For  the  GS  theory  to  be  validated,  accurate  and  meaningful  estimates  of  the  three  GS  parameters  must  be  achieved.  It  is  not 
likely  that  this  can  be  accomplished  with  field  data.  Laboratory  conditions  will  likely  provide  better  control  over  sediment 
properties  and  more  accurate  measurements.  But  this  analysis  of  a  large  set  of  field  observations,  coupled  with  the  sensitivity 
analysis,  provides  a  first  step  in  bounding  the  range  of  values  that  should  be  investigated. 

The  estimate  of  w  is  a  critical  first  step  and  relies  on  accurate  measurement  of  porosity,  compressional  sound  speed  and 
attenuation,  all  of  which  contribute  significantly  to  the  error.  It  is  clear  from  this  data  set  that  a  single  value  of  n  cannot  be 
assumed  for  all  sediments.  But  more  accurate  measurements  are  required  in  order  to  determine  whether  these  GS  parameters  can 
be  specified  so  that  GS  theory  can  be  thoroughly  validated. 

APPENDIX  -  PARTIAL  DERIVATIVES 


List  of  symbols 

Co  compressional  sound  speed  in  equivalent  suspension 
Cp  compressional  wave  speed  in  sediment 

Cs  shear  wave  speed  in  sediment 

Kq  Bulk  modulus  of  equivalent  suspension 

Bulk  modulus  of  grain 
Bulk  modulus  of  pore  fluid 
n  Strain-hardening  index 

N  Porosity 

T  Arbitrary  time  constant 

OCp  Compressional  attenuation  coefficient 

^  Dimensionless  grain-shearing  coefficient 

Yp  Compressional  wave  grain-shearing  modulus 

Ys  Shear  wave  grain-shearing  modulus 

Po  Bulk  density  of  sediment 

(O  Angular  Frequency 

The  partial  derivatives  of  the  Buckingham  model  parameters  with  respect  to  measured  values  are  given  below.  For  the  sake  of 
clarity,  some  derived  parameters  are  given  first.  Then  relevant  expressions  are  given  in  terms  of  those  derived  parameters. 


|L= ^  si„(„*)+2X  si„=  (<? )  I  - 

sin(.:r)^-2Xsln^(‘f)H-2^)‘(ux‘) 

— =  — (^tsinV^Yl  +  X^f  -^--^ 

^  ;cvV-/  ''2  A  y  ^ 


TtK 


sin(n;r)+ 2X  sin^  1  +  2 


^=-4lt.g(s«')+fcot(«)]+2 


3h 


(n-x‘) 


%  =  -f  wJW  si"(“)-  '"SW 


_  2 

^  <?te 


=  ^K)''  t  +  |i+^0fiF)''J^(log(ft;7)+f  vy 


^=^(»r)"  s  l.xiJc^-fr  §.2-3  |i.2-(y<»r)-f  (,og(»r).f;>- 


<?<Qr 

_ j 

2cq 

|^=^K)"S  [+2'0<»r)"f  /  ^ 


2poCo  ^  ^ 


Acknowledgment 


The  authors  would  like  to  thank  Michael  Buckingham  for  sharing  the  MATLAB  code  for  his  GS  model.  This  work  was 
supported  in  part  by  the  Office  of  Naval  Research. 


References 

[1]  M.  J.  Buckingham,  “Wave  propagation,  stress  relaxation,  and  grain-to-grain  shearing  in  saturated,  unconsolidated  marine  sediments,” 7.  AcousL  Soc.  Am.^  108 

(6),  pp.  2796-2815,  2000. 

[2]  M.  A.  Biot,  “Theory  of  propagation  of  elastic  waves  in  a  fluid-saturated  porous  solid  I.  Low  frequency  range,”  J.  Acoust.  Soc.  Am..,  28  (6),  pp.  168-178,  1956. 

[3]  M.  J.  Buckingham,  “On  pore-fluid  viscosity  and  the  wave  properties  of  saturated  granular  materials  including  marine  sediments,”  7.  Acoust.  Soc.  Am.,  122  (3), 

pp.  1486-1501,2007. 

[4]  D.  R.  Jackson  and  M.  D.  Richardson,  High-Frequency  Seafloor  Acoustics,  Springer,  pp.  123-170,  2007. 

[5]  A.  B.  Wood,  A  Textbook  of  Sound,  third  ed.,  G.  Bell,  1964. 

[6]  M.  J.  Buckingham,  “Compressional  and  shear  wave  properties  of  marine  sediments:  Comparisons  between  theory  and  data,”  J.  Acoust.  Soc.  Am.,  1 17  (1),  pp. 

137-152,  2005. 

[7]  B.  T.  Heffner  and  K.  L.  Williams,  “Sound  speed  and  attenuation  measurements  in  unconsolidated  glass-bead  sediments  saturated  with  viscous  pore  fluids,”/ 

Acoust.  Soc.  Am.,  120  (5),  pp.  2538-2549,  2006. 

[8]  A.  Tarantola,  Inverse  Problem  Theory:  Methods  for  Data  Fitting  and  Parameter  Estimation,  Elsevier,  pp.  68-69,  1987. 

[9]  E.  I.  Thorsos,  K.  L.  Williams,  N.  P.  Chotiros,  J.  Christoff,  K.  Commander,  C.  F.  Greenlaw,  D.  V.  Holliday,  D.  R.  Jackson,  J.  L.  Lopes,  D.  E.  McGehee,  J.  E. 

Piper,  M.  D.  Richardson,  and  D.  Tang,  “An  overview  of  SAX99:  Acoustic  measurements”,  lEEEJ.  Ocean.  Eng.,  26,  pp.  26-53,  2001. 


